

//#include <feel/options.hpp>
//#include <feel/feelcore/application.hpp>

#include <feel/feelfilters/geotool.hpp>



int main(int argc, char** argv)
{
    using namespace Feel;

    double meshSize=0.05;

#if 0
    //Geometry
    GeoTool::Node x1(0,0);
    GeoTool::Node x2(4,1);
    GeoTool::Rectangle R( meshSize,"LEFT",x1,x2);
    R.setMarker(_type="line",_name="Entree",_marker4=true);
    //R.setMarker(_type="line",_name="Paroi",_marker1=true,_marker3=true);
    R.setMarker(_type="line",_name="ParoiB",_marker1=true);
    R.setMarker(_type="line",_name="ParoiH",_marker3=true);
    R.setMarker(_type="line",_name="Sortie",_marker2=true);
    R.setMarker(_type="surface",_name="OmegaFluide",_markerAll=true);


    typedef Simplex<2,1,2> convex_type;
    typedef Mesh<convex_type> mesh_type;

    auto mesh = R.createMesh<mesh_type>("geotest1d");
#else

    //Geometry
    GeoTool::Node x1(0,0.5);
    GeoTool::Node x2(5,0.5);
    GeoTool::Line R( meshSize,"LEFT",x1,x2);
    /*R.setMarker(_type="line",_name="Entree",_marker4=true);
    //R.setMarker(_type="line",_name="Paroi",_marker1=true,_marker3=true);
    R.setMarker(_type="line",_name="ParoiB",_marker1=true);
    R.setMarker(_type="line",_name="ParoiH",_marker3=true);
    R.setMarker(_type="line",_name="Sortie",_marker2=true);
    R.setMarker(_type="surface",_name="OmegaFluide",_markerAll=true);*/


    typedef Simplex<1,1,2> convex_type;
    typedef Mesh<convex_type> mesh_type;

    auto mesh = R.createMesh<mesh_type>("geotest1d");
#endif

    return 0;

}


